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The alternative dynamics of loop quantum cosmology is examined by the path integral formulation. We 
consider the spatially fiat FRW models with a massless scalar field, where the alternative quantization inherit 
more features from full loop quantum gravity. The path integrals can be formulated in both timeless and depa- 
rameterized frameworks. It turns out that the effective Hamiltonians derived from the two different viewpoints 
are equivalent to each other. Moreover, the first-order modified Friedmann equations are derived and predict 
quantum bounces for contracting universe, which coincide with those obtained in canonical theory. 
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O ! I. INTRODUCTION 

(N 

In recent 25 years, considerable progress has been made in loop quantum gravity (LQG)[l-4], which is a non-perturbative 
and background independent quantization of general relativity. A particular feature of LQG is to to inherit the core concept of 
general relativity that the fundamental physical theory should be background independent. The construction of LQG inspired 
the research on spin foam models (SFMs) [5], which are proposed as the path-integral formulation for LQG. In SFMs, the 
transition amplitude between physical quantum states is formulated as a sum over histories of all physically appropriate states. 
The heuristic picture is the following. One considers a 4-dimensional spacetime region M bounded by two spacelike 3-surfaces, 
which would correspond to the initial and final states for the path integral. The feature of SFMs is to employ the spin network 
' states in LQG to represent the initial and final states. A certain path between the two boundaries is a quantum 4-geometry 
interpolated between the two spin networks. The interpolated spin foams can be constructed by considering the dual triangulation 
of M and coloring its surfaces with half integers and edges with suitable intertwiners. In order to obtain the physical inner product 
i—{ ' between the two boundary states, one has to sum over all possible triangulations and colorings [1]. This is certainly difficult to 
J> , be achieved. Thus it is desirable to examine the idea and construction of SFMs by certain simplified models. 

It is well known that the idea and technique of canonical LQG are successfully carried out in the symmetry -reduced models, 
known as loop quantum cosmology (LQC) [6]. Because of the homogeneity (and isotropy), the infinite degrees of freedom of 
■ gravity have been reduced to finite ones in LQC [7]. Hence LQC provides a simple arena to test the ideas and constructions of the 
full LQG. Therefore, it was proposed to test the idea and construction of SFMs by LQC models [8]. It was shown in Refs.[9, 10] 
that the transition amplitude in the deparameterized LQC of k = Friedmann universe equals the physical inner product in 
the timeless framework, and concrete evidence has been provided in support of the general paradigm underlying SFMs through 
the lens of LQC. How to achieve local spinfoam expansion in LQC was also studied in Refs. [11, 12]. Recently, the effective 
Hamiltonian constraint of k = LQC was derived from the path integral formulation in timeless framework in Ref.[13], and the 
effective measure introduces some conceptual subtleties in arriving at the WKB approximation. The path integrals and effective 
Hamiltonians of LQC for k = +1,-1 FRW models are also derived [14]. The coherent state functional integrals are also 
proposed to study quantum cosmological models [15]. However, in canonical LQC, there are some quantization ambiguities in 
• constructing the Hamiltonian constraint operator. Theoretically, the quantization method which inherits more features from full 
LQG would be preferred. In spatially fiat FRW model, two alternative versions of Hamiltonian constraint operator are proposed 
in Ref. [16], where the big bang singularity resolution, correct classical limit and re-collapse of an expanding universe because 
of higher-order quantum effect are obtained. 

The purpose of this paper is to study the alternative effective dynamics of LQC from the viewpoint of path integral. The 
models which we are considering is the spatially flat FRW models with a massless scalar field. We will formulate the path 
integrals for two different proposed dynamics of LQC in both timeless and deparameterized frameworks. The multiple group 
averaging method proposed in Ref. [ 14] is employed. It turns out that the alternative effective Hamiltonians derived from the two 
different viewpoints are equivalent to each other. Moreover, the first-order modified Friedmann equations can be derived, and 
the quantum bounces for contracting universe will also be obtained in both kinds of dynamics. In section 2, we will introduce 
the basic frameworks of LQC and the path integral approach of multiple group averaging. Then we will derive the path integral 
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formulation and effective Hamiltonian of alternative quantization model I and model II in sections 3 and 4 respectively. Finally 
a summary will be given in section 5. 

n. BASIC SCHEME 

We are consider the spatially fiat FRW universe filled with a massless scalar field. In the kinematical setting, it is convenient to 
introduce an elementary cell V in the spatial manifold and restricts all integrations to this cell. Fixing a fiducial fiat metric °q ab and 
denoting by V the volume of V in this geometry, the physical volume reads V — a 3 V . The gravitational phase space variables 

— (1/3^ — (2/3) 

— the connections and the density weighted triads — can be expressed as = cV °u l a and Ef = pV a \f% °ef, 
where {°(jJ l a , °ef) are a set of orthonormal co-triads and triads compatible with °q ab and adapted to V. p is related to the scale 

factor a via \p\ = Vo^a 2 . The fundamental Poisson bracket is given by: {c, p} = SirGj/S, where G is the Newton's constant 
and 7 the Barbero-Immirzi parameter. The gravitational part of the Hamiltonian constraint reads C grav = — &c 2 y/\p\/^ 2 . To 
apply the area gap A of full LQG to LQC [17], it is convenient to introduce variable p, satisfying 

p 2 \p\ = A = (WSrry)^, (2.1) 

where £ 2 = Gh, and new conjugate variables [16, 18-20]: 

;= s g n(p)\p\i b ^pc 
2n 1 e 2 VA ' 2 

The new canonical pair satisfies {b, v} = 4. On the other hand, the matter phase space consists of canonical variables <f) and 
which satisfy {4>,p$} = 1. To mimic LQG, the polymer-like representation is employed to quantize the gravity sector. 
The kinematical Hilbert space for gravity then reads HE™ = £ 2 (RBohr, dpH ), where Reohr is the Bhor compactification of 
the real line and dpn is the Haar measure on it [21]. It turns out that the eigenstates of volume operator v, which are labeled 
by real number v, constitute an orthonomal basis in "Hf[^ v as (vi\v-2) = 5 Vl:V2 . For the scalar matter sector, one just uses the 
standard Schrodinger representation for its quantization, where the kinematical Hilbert space is "HJJJn" = L 2 (M., d<p). The total 
kinematical Hilbert Space of the system is a tensor product, %km = ^kirT ® ^kiSf*' °f tne arjove two Hilbert spaces. As a 
totally constrained system, the dynamics of this model is reflected in the Hamiltonian constraint C grav + Cmatt = 0. Quantum 
mechanically, physical states are those satisfying quantum constraint equation 

(Cgrav + Cmatt)*(«, 4>) = 0, (2.3) 

which is not difficult to be rewritten as an Klein-Gordon like equation [22]: 

C*(«,0) = (^-e)*(i; I ^)=O. (2.4) 
Eq.(2.4) indicates that we can get physical states by group averaging kinematical states as [11] 

da e lad f(v, 0), V/(«, 4>) G H ki „, (2.5) 

and thus the physical inner product of two states reads 

( /| 5 ) phy = (V f \g ) = lim r da (f\e lad \g). (2.6) 
As is known, in timeless framework the transition amplitude equals to the physical inner product [9, 10], i.e., 

/Oio _ 

da(vf,4>f\e laC \vi 7 (pi). (2.7) 
-OCo 

On the other hand, Eq.(2.4) can also be written as 

fl$tf(w, 4>) + <f) = 0. (2.8) 
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The similarity between Eq.(2.8) and Klein-Gorden equation suggests that one can regard 4> as internal time, with respect to which 
gravitational field evolves. In this deparameterized framework, we focus on positive frequency solutions, i.e., those satisfying 



-id^+(v, <t>) = ve*+(«, 0) = m+(v, 0). (2.9) 

The transition amplitude in deparameterized framework is then given by 

A dep (v f ^ f ; v h ct>i) = (Vfl^f-^lvi), (2.10) 

where \vi) and \v/) are eigenstates of volume operator in 'Hfj„ v , and <fi is the internal time. 

Starting with (2.7), we now consdier the transition amplitude of path integral under timeless framework. In order to compute 
(vf, 4>f\e laC \vi, <pi), a straightforward way is to split the exponential into N identical pieces and insert complete basis as in [13]. 
However, since a is the group-averaging parameter which goes from — oo to oo, it is unclear whether a could be treated as the 
time variable t in non-relativistic quantum mechanics path integral. We thus consider alternative path integral formulation of 
multiple group-averaging [14]. Here one single group averaging (2.5) is generalized as multiple ones: 

lim / dae la5 \v,(l)) 



I™ 77^ — / d&N...—^ d& 2 e l( " 1+ - +<ijv)c |t>,0). (2.11) 



a No ,...,a lo ^oo 2a No J_ 5lNo 2a 2o J-a 2o J - aio 

In order to trace the power for expansion, we re-scale the parameters by a n = ea n , where e = -h. Then (2.7) becomes 

Au s {vf,(j)f) Vi,(j)i) 

1 r aN - i r a *> r a >° % y. 

= lim / da N ...- / da 2 ■ e / da>i(vf, (j)r\e " =1 \vi,<f>i), (2.12) 

a N o,-,<xi„->oc 2a No J_ aelo 2a 2o J- a2o J- a , a 

where a„ = a„„/e, n — 1,2, .., N. Next, we are going to insert a set of complete basis at each knot. Notice that \v, <fi) is the 
eigenstate of both volume operator and scalar operator simultaneously in Hkm, which is written as \v)\(j)) for short, and 



(2.13) 



Thus, we have 

N 



(v f ,4>f\e£™ nC \v i ,4>i)= J2 [ d(f>N-i...dct> 1 f[((t> n \{v n \e iea « d \v n - 1 )ct> n _ 1 ), (2.14) 

UjV_ i , . . .vi n—1 

where Vf = Vn, <\>f = 4>Ni v% — vq, <f>i = (f)Q have been set. Since the constraint operator C has been separated into gravitational 
part and material part, which live in %j^ v and separately, we could calculate the exponential on each kinematical space 

separately. For the material part, one gets 

■ £| f iea f\ 

(4> n \e tan ~iP\4> n ^i) = / dp 4>n {(t) n \p rt>n ){p 4>n \e' ian '^\(t)n-i) 

W**- e T ^ ^' <2j5) 

As for the gravitational part, in the limit N — >• oo(e — >• 0), the operator e ^ lea « e C an be expanded to the first order, and hence we 
get 

(v n \e-' lea ^\v n ^) = ^„,„„_ 1 - iea n {v n \Q\v n ^) + 0(e 2 ). (2.16) 



III. PATH INTEGRAL AND EFFECTIVE DYNAMICS OF MODEL I 



In this section, we employ one of the alternative Hamiltonian constraint operators for LQC proposed in [16], where the 
extrinsic curvature K l a in the Lorentz term of the gravitational Hamiltonian constraint was quantized directly following the 
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procedure in full LQG. Using the simplification treatment in [18], we can get the action of the gravitational Hamiltonian operator 
in this model as: 

3ttG 7 



O f \v) = ^^v[(v + 2)\v + 4) - 2v\v) + (v-2)\v- 4)] 



3ttG(1 + 7 2 ) 
16 



v[(v + 4)|u + 8) - 2v\v) +(v- A)\v - 8)] . 



(3.1) 



which leads to 



(« R |e F K_i) 



3nGj 2 v„ + v n . 



-v n -i- 



— {3v n ,v n -i+4 ~ 23v n ,v n -i + $v n , v n -i— 4) 



37rG(l +7 2 ) Vn + Vn-i , x ox , x , 

Vn-1 ^ \°v n ,V n -l+8 ~ ZO i>n^n-l "+" "i>„,u„_i-8j. 



(3.2) 



Applying (3.2) to (2.16) and writing Kronecker delta as integral of b n , which acts as the role of conjugate variable of v n , we 
have 



(v n \e- Ua " e °\v n ^) 



2 n 



db n e 



-«„(i>„-«.-i) fi_ ia^S^Vn 



V n + V n -1 . 2 



sin 2 (26 n )[l - (1 + 7 2 ) sin 2 (26 n )]) + 0(e 2 ). 



(3.3) 



Applying (2.15) and (3.3) to (2.12), and then taking 'continuum limit', we obtain 



lim 

N—>oo a N , 



- ~ ( e n 2a „ o )/ 



do 



N ■•■ 



— a No J —ctio 



da 1 



x d<f> 



(—) 

\2Tih 



N 



dp < f >N ...dp<i> 1 ^2 (-) / db N ...dbi 



N 

x J^J exp ie 

n=l 



VN-1, — ,V1 

/ m 2 



N 



P<f> n <Pn ~ <Pn-l , Vn ~ «m-l , 

-r o« h a n -hr - i>nGv n -\ - 

he e \ fi z 2 



^■sin 2 (26„)[l-(l+ 7 2 )sm 2 (2b n )] 



(3.4) 



Finally, we could write the above equation in path integral formulation as 



A *\s{Vf,<t>f, Vi,(f>i) 



--c I Da I V(f> I V P4> I Vv I Vb cxp - / dr 



Pit -hbv + haljjt- 3ttGv 2 sin 2 (26) [1 - (1 + 7 2 ) sin 2 (26)] 



(3.5) 



where c is certain constant, and a dot over a letter denotes the derivative with respect to a fictitious time variable r. The effective 
Hamiltonian constraint can be read out from Eq. (3.5) as 



C F ff = p - SttGv 2 sin 2 (26) [1 - (1 + 7 2 ) sin 2 (26)] 



(3.6) 



Using this effective Hamiltonian constraint, we can explore the effective dynamics of universe by the modified Friedmann 
equation: 



m = 



SirGp 7 2 + 4(1 + 7 2 WPc 7 2 Pq f, 4(1 + 7 > 

3 I I + 7 2 + 2(l + 7 2 ) 2 p 



1 - a/1 



4(1 + 7 > 



Pc 



(3.7) 



where p = -J^? is the matter density and p c = 327r glft 7 :i is a constant. This modified Friedmann equation coincides with the 
one in [16] if we ignore the higher-order quantum corrections therein. It is easy to see that if the matter density increase to 
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P = 4(x+ 7 5) ' me Hubble parameter would be zero and the bounce could occur for a contracting universe. On the other hand, in 
the classical region of large scale, we have p <C p c and hence Eq. (3.7) reduces to the standard classical Friedmann equation: 

tt2 _ 8ttGp 
n 3 . 

Besides the timeless framework, we can also employ the above group-averaging viewpoint for the deparameterized framework. 
To this end, we define a new constraint operator C+ = t£ — %. Then Eq.(2.9) can be rewritten as 

Cly + (v,<p)=Q. (3.8) 
The transition amplitude for this new Hamiltonian constraint reads 

da{v J ,tf>f\e iaC +\v i ,tf>i)= lim / da(v f , (/)f\2\p 4> \8(p 4> )e iaC \v l , fa), (3.9) 

-a„ a °^ roc J-a 

where 

WK) = KIK> ; ^)K> = ||^ } (3-io) 

Similar to the timeless case, the integration over single a can be written as multiple integrations, 

Aiep(v/, <t>f, Vi,<pi) 



1 f aN ° 1 [ a2 ° f a '° i £ ea„C 



= lim / da N ...- / da 2 ■ e dai (v f , <fi f \e ™ =1 </>»)■ (3.11) 

a„„,...,Q,„^oo 2a No J_ am 2a 2o J- a2o J- ai0 

Since we work in the deparameterized case, it is reasonable to insert in the completeness relation [23] 
where 

|v, </>,+)= lim / dP e tf>c +\v,<t>). (3.13) 

Note that the Hilbert space in deparameterized framework is unitarily equivalent to the physical Hilbert space with the basis 
(3.13). Then the first piece of the exponential in (3.11) becomes 

/ceio 
d(eai) (ui,0i,+|e ieaiC +|uo,0o) 
-aio 

rP'io , - 

= lim / dft (vi,<t>! I e ift C+ 1 fo^o) 

= lim e / dA<ui,0i|e ,e * o +|w(,,0o}, (3.14) 

and the last piece of the exponential reads 

lim / da N / d0' N (v I r,<l> N \e ieai ' c +e i l i » c +\v f r-i,<f> N -i) 

a No ,P' No ^oo Za No J^ aNo J-l3' No 

= lim / dj3' N (ujv,</>7v|e l/3 " c +|«Ar-i,0w-i) 

= lim e f dPNivN^Nle^^Vif-uto-i). (3.15) 
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The remaining pieces of the exponential can also be expressed as 

1 



lim 

a„„- >oo 2a 



lim 



= lim e 



no J -a 
P'n 



da n {v n ,(f> n , +| e 8ea ™ c '+ \v n -x, n _i,+) 



|8„C 



d/3 n (w n ,^ n |e^™ c +|f„_i,0 n _i). 



(3.16) 



So (3.11) becomes 

lim e w 

f3 No ,...,()i -+oo 



dj8 



JV • 



JV 



-0i 



dfii [ K>n|e^" c+ k-i,0„-i) 



= lim 



JV 



(9« 



d/3 



JV 



d/3i [ K,0n|2b J0(p Je le ' 3 " c | V „_i,^-i), (3.17) 



•UjV — 1 ,■ ■ ■ ,^1 Tl— 1 



where (3.9) is applied in second step. Now we can split each piece in (3.17) into gravitational and material parts. Calculations 
similar to those in timeless framework lead to 



^dep(>/,0/; Vi,(j>i) 



JV 



: lim lim 

JV— >oo /3no,'",Pio— >"00 
JV 

x JJ 2|p ?iri |6'(p0 ll )expie 



d/3 



JV- 



d/3; 



\2llh) 



JV 



dp rt>N ...dp <t>1 ^2 



N 



dbisr...dbi 



v n - u„_i 



ft e 

,,2 



+ /Jn(^ - 3^ V" ^ 1 sin 2 (26 n )[l - (1 + 7 2 ) sin 2 (26„)]) 
We can integrate out /3„ and and arrive at 



(3.18) 



^dep (Vf,<t>f, V l ,(j) l ) 

lim y 



JV 



DJV_1,...,«1 



JV 



db]y...dbi 



JV 



I expie(cf>f 



V n + V n -l . 2 



sin 2 (26„)[1 - (1 + 7 2 ) sin 2 (26„)] - bj n 



'e(0/ - 00 



c' / Vv I Vb exp| [ #(y 'z-kGW-v 1 sin 2 (26) [1 - (1 + 7 2 ) sin 2 (26)] - fiufc), 



where v = ^ . The effective Hamiltonian in deparameterized framework can be read out from (3.19) as: 



ft^ = -^3nGK 2 v 2 sin 2 (26) [1 - (1 + 7 2 ) sin 2 (26)], 
We can use this effective Hamiltonian to get the effective equations of motion: 

SirGv 2 sin (46) [1 - 2(1 + 7 2 ) sin 2 (26)] 



yj 3irGv 2 sin 2 (26) [1 - (1 + 7 2 ) sin 2 (26)] 



6 = -sgn(v) J ZirGv 2 sin 2 (26) [1 - (1 + 7 2 ) sin 2 (26)], 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



which also predict a bounce of a contracting universe when v = 0. It is easy to see that the bounce point coincides with that from 
(3.7). In fact, the effective equations derived in timeless and deparameterized frameworks coincides with each other. Hence, at 
least at first-order level, both methods of path integral confirm the effective dynamics of this model in canonical theory. 
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IV. PATH INTEGRAL AND EFFECTIVE DYNAMICS OF MODEL II 



In the other Hamiltonian constraint operator proposed in [16], the Lorentz term was constructed by using the fact that the 
extrinsic curvature K l a is related to the connection A l a by A % a = 7 A"* in the spatially fiat case. Reexpressing the connection by 
holonomy, we can get a simplified version of this operator (9 s similar to F as: 

37rG 7 2 



-v[{v + 2)\v + 4) - 2v\v) + (v - 2)\v - 4)] 



- 3ttG(1 + j 2 )v[(v + l)\v + 2) - 2v\v) + (v - l)\v - 2)] . 
Following the same procedure in the previous section, we can get the transition amplitude in timeless framework as: 

A t\ s ( v f^f, v %Ai) 



(4.1) 



lim 

/V->oc 



lim 

CKA'o,...,CK/ — ^OO 



N 1 



(—) 

V 27T?l/ 



N 



da 



N '•■ 



da.\ 



E 



N 



^ cxp t 



P0„ 0n - 



V n -1 



which gives an effective Hamiltonian constraint: 



C c s ff = S - 127rGi; 2 sin 2 6(1 + 7 2 sin 2 6) 
From this effective constraint, we can derive another modified Friedmann equation: 

3(l+ 7 2 ) + 7 2 p/Pc , 2(l + 7 2 )Pc ' 



vn-i,---,-vi 



db]y...dbi 



l2nGv n 



V n -l 



sin 2 6„(1 + 7 2 sin 2 6„) 



Hi = 



8nGp 



1 - 



r 



7 4 p 



1 



7 2 P 



1 



2 

Pc 



- 1 



(4.2) 



(4.3) 



(4.4) 



It is obvious that this effective equation is different from Eq.(3.7) of the model I. By Eq.(4.4), the quantum bounce would occur 
when matter density increase to p = 4(1 + "f 2 )p c for a contracting universe. It is easy to see that the Friedmann equation (4.4) 
can also reduce to the classical one when p <C p c - Moreover, Eq.(4.4) coincides with the corresponding effective Friedmann 
equation in [16] if the higher-order quantum corrections therein are neglected. 
Similarly, we can also get the deparameterized amplitude: 



lim 



N- 



l)JV_l,...,«l 



dbN---db\ 



N 



J expie(0 / - fa) 



71=1 



-1 ~ sm b n {l + Y sm b n )-b n — 



\2ixGv n 



v n -i 



e((f>f - 00 



=d J Vv J Vb ex P^y d< t>{\J 12nGh 2 v 2 sin 2 6(1 + 7 2 sin 2 b) - hvb), 
which gives an effective Hamiltonian: 

Hi 



- v 'l2TrGh 2 v 2 sin 2 6(1 + 7 2 sin 2 6). 

The evolution of v and 6 with respect to <fr can be obtained from this effective Hamiltonian as: 

. _ dv _ \2ixGv 2 sin (26) (1 + 2 7 2 sin 2 6) 
V= ~d4>~ 



(4.5) 



(4.6) 



2J12ttGH 2 v 2 sin 2 6(1 + 7 2 sin 2 6) 



6 = 



db 



sgn(w) J \2vGh 2 v 2 sin 2 6(1 + 7 2 sin 2 6), 



(4.7) 



which again coincide with the effective Eq.(4.4) and predict the same bounce when v = 0. Hence, the corresponding effective 
dynamics of this model in canonical theory is also confirmed at first order by the path integral. 
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V. SUMMARY 

The motivations to study alternative dynamics of LQC are in two folds. First, since there are quantization ambiguities in 
constructing the Hamiltonian constraint operator in LQC, it is crucial to check whether the key features of LQC, such as the 
quantum bounce and effective scenario, are robust against the ambiguities. Second, since LQC serves as a simple arena to test 
ideas and constructions induced in the full LQG, it is important to implement those treatments from the full theory to LQC as 
more as possible. Unlike the usual treatment in spatially flat and homogeneous models, the Euclidean and Lorentz terms have 
to be quantized separately in full LQG. Therefore, this kind of quantization procedure which kept the distinction of the Lorentz 
and Euclidean terms was proposed as alternative dynamics for LQC [16]. It was shown in the resulted canonical effective theory 
that the classical big bang is again replaced by a quantum bounce. Hence it is desirable to study such kind of prediction from 
different perspective. Moreover, it is desirable to examine the idea and construction of SFMs by the simplified models of LQC. 

The main results of the present paper can be summarized as follows. The path integral formulation is constructed for spatially 
flat FRW models under the framework of LQC with two alternative dynamics. In both models, we can express the transition 
amplitude in both the timeless and the deparameterized frameworks by multiple group averaging procedure. We can derive the 
effective Hamiltonians from both viewpoints. It turns out that in both models the resulted effective dynamics from the timeless 
and the deparameterized path integrals are equivalent to each other. This indicates the equivalence between two frameworks 
of path integral. Moreover, the modified Friedmann equations for both models are also obtained and coincide with the cor- 
responding equations in [16] if the higher-order quantum corrections therein are neglected. This indicates the equivalence of 
the canonical approach and the path integral approach in LQC. Since our path integral approach inherits significant features 
of SFMs, it provides certain good evidence to support the scheme of SFMs. In both models, the quantum bounce will replace 
the big bang singularity due to the modified Friedmann equations when the matter density increase to the magnitude of Planck 
density. Hence the quantum bounce resolution of big bang singularity in LQC is robust against the quantization ambiguities 
of the Hamiltonian constraint. Moreover, the alternative modified Friedmann equations (3.7) and (4.4) set up new arenas for 
studying phenomenological issues of LQC. 
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